Scintillation in the Circinus Galaxy H2O megamasers 
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ABSTRACT 

We present observations of the 22 GHz water vapor megamasers in the Circinus galaxy made 
with the Tidbinbilla 70 m telescope. These observations confirm the rapid variability seen earlier 
by Greenhill et al. (1997). We show that this rapid variability can be explained by interstellar 
scintillation, based on what is now known of the interstellar scintillation seen in a significant 
number of flat spectrum AGN. The observed variability cannot be fully described by a simple 
model of either weak or diffractive scintillation. 

Subject headings: galaxies,: active, - galaxies: individual (Circinus) - galaxies: Seyfert - ISM: molecules 
- masers: ISM: interstellar scintillation 



1. Introduction 

Water vapor megamaser emission was first de- 
tected from the Circinus galaxy by Gardner & 
Whiteoak (1982), making it the second closest 
Seyfert galaxy, after NGC 4945, to host water va- 
por megamaser emission. Recent VLBI observa- 
tions show the emission arises in a thin warped 
disc with a clumpy, wide-angle outflow subtending 
an angle of ~ 80 mas which corresponds to a lin- 
ear region of ~ 2 pc at a distance of 4.2 Mpc. The 
nucleus of the galaxy hosts a 1.7 ± 0.3 xlO^ Mq 
black hole which is also a powerful X-ray source 
(GreenhiU et al. 2003). 

Greenhill et al. (1997) observed very rapid 
variations in several of the strong maser lines, 

^Affiliated with the Australia Telescope National Facil- 
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with timescales as short as a few minutes in one 
line. They discussed both intrinsic variability and 
diffractive interstellar scintillation (DISS) as po- 
tential mechanisms, and concluded that both were 
possible although favouring DISS. Both mecha- 
nisms placed limits on the physical conditions 
of the masers, with intrinsic variations requiring 
short gain paths of 1-10 AU, and DISS requiring 
angular sizes < 1.75 //as (< 7 AU at a distance of 
4.2 Mpc). The implied brightness temperatures 
are in excess of 10^^ K. Similar conclusions were 
reached by Maloney (2002). 

Over the last five years unequivocal evidence 
has accumulated that weak interstellar scintilla- 
tion (WISS) is the principal mechanism responsi- 
ble for the intra-day variability (IDV) observed at 
cm-wavelengths in many compact, flat-spectrum 
AGN (e.g Jauncey et al.(2003) and references 
therein) . In the light of these recent developments 



we show here that the very rapid variability seen 
in Circinus is also well described by ISS. 

2. Observations 

We present data from three epochs, one obser- 
vation made with the Parkes 64 m telescope on 
DOY 289 of 1995 and two observations made in 
1996 with the Tidbinbilla 70 m. The Parkes obser- 
vation has been previously analysed and published 
by Greenhill et al.(1997). The details of the ob- 
servation and data reduction are not reproduced 
here. 

The Tidbinbilla 70 m observations were made 
at high elevations in left circular polarization on 
1996 March 4 and March 28 (DOYs 124 and 
148) using a 20 MHz bandpass with 1024 spectral 
channels giving a channel spacing of 19.5 KHz or 
0.26 kms"-*^. The bandpass was centred on 537 
kms~^with a velocity range of 402-672 kms^^on 
DOY 124, and centred on 355 kms~^with a ve- 
locity range of 220-490 kms'^on DOY 148 (Helio- 
centric frame of reference). Circinus was observed 
in position switching mode with groups of eight 
45-second integrations on- and off-source. The to- 
tal duration of the observations was 2.8 hours on 
DOY 124 and 1.9 hours on DOY 148. 

The data on both days were calibrated assum- 
ing a constant 150 Jy system temperature so there 
may be small, slow trends due to atmospheric ef- 
fects and slow changes in gain with elevation. The 
weather on both occasions was fine and the mea- 
sured gain-elevation curve shows that such effects 
will be both slow and small. Also, there is no 
significant correlation between the observed vari- 
ations of different maser features which indicates 
that our calibration is adequate. The measured 
RMS noise level in a spectral channel in a 45 sec- 
ond integration was 0.11 Jy and 0.17 Jy on DOY 
124 and 148 respectively. The average spectrum 
is shown in Figure 1. 

Characteristic time-scales for variability, Tchar, 
were calculated from the normalized auto-correlation 
function (ACF) of the light curves for both the 
Parkes and Tidbinbilla data, using the discrete au- 
tocorrelation function (Edclson & Krolik 1988). 
Following Rickett et al.(2002), we measure the 
lag at which the ACF falls to 0.5 as the estimate 
of Tchar- The velocities of the Parkes data have 
been re-calculated as it appears that these were 



previously in reference to the LSR despite being 
labeled as heliocentric in Greenhill et al. (1997). 

There are some difficulties in calculating the 
ACF of the data, due to the regular observing 
gaps. This is because the number of points for 
a lag interval decreases as the lag time approaches 
the gap time. The regular sampling pattern, with 
equal integration times on- and off-source and 
losses due to calibration and slewing overheads 
means that there are no data for some lags. Using 
the discrete correlation function, we have excluded 
those lags that were computed using less than a 
third of the available data. 

3. Results and Discussion 

Table 1 summarises the measured properties of 
the well-observed lines in the data. The velocity 
and FWHM of the lines arc taken from the fit- 
ted Gaussian profiles as is the mean flux density 
(< S >) and standard deviation of the variariabil- 
ity (ct). The modulation index, /i, was calculated 
by cr/ < S >. The characteristic timescale, as 
defined above, is included together with the es- 
timated depth of the ACF at its first minimum. 
The listed errors were calculated using the esti- 
mation error formula of Appendix A of Rickett 
et al. (2002). The extremely large errors are due 
to the short observing time and illustrate the need 
for longer observations. 

The light curves presented in Figures 2 and 3, 
as well as the original Parkes light curves (Green- 
hill et al. 1997), clearly show the rapid variability 
in the stronger lines. It is apparent that there is 
no significant cross-correlation in the Tidbinbilla 
data between the variations at the different line 
frequencies, showing that the variations cannot be 
due to systematic effects, miss-pointing or atmo- 
spheric effects for example. Moreover, different 
lines show quite different time scales, with the 295 
kms~^line on DOY 148 showing much slower vari- 
ations than the 257 kms~^line. There is, however, 
some evidence for rapid, low amplitude variations 
superimposed on the slow variations in the 295 
kms~^line. 

There are two strong lines present in both the 
Parkes and Tidbinbilla observations (the 295 and 
557 kms~^lines). In both cases, they are much 
weaker in the Tidbinbilla observations, having 
roughly halved their flux densities. The modu- 
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lation indices, timcscales and FWHM (Full-width 
half-maximum of the fitted Gaussian profile) of 
the masers have also changed. In the DOY 124 
observations of the 557 kms~'^line, it has narrowed 
considerably and has a slightly faster timescale 
while the modulation index is very similar. For the 
295 kms~^line, its FWHM is unchanged while the 
modulation index has doubled and the timescale 
has significantly increased. A prediction of ISS 
theory is that the characteristic timescale of vari- 
ability should change throughout the year as the 
motion of the Earth around the Sun changes the 
relative speed of the ISM, across the line of sight 
(Macquart & Jauncey 2002). This should then 
cause the variability of both lines to speed up 
or slow down. The observed timescale changes, 
where one line slows down dramatically and the 
other speeds up tells us that source structure is 
affecting the scintillation and may be evolving on 
timescales of a few months. The observed line nar- 
rowing in the 557 kms~^line is indicative of this 
and previous observations have indicated that the 
average lifetime of an individual maser feature is 
^ 1 year. 

3.1. Weak scintillation 

The quasi-sinusoidal variations of the light 
curves in Figures 2 and 3 plus those in Green- 
hill et al's Figure 2, are remarkably similar in 
amplitude and time scale to the variations seen 
at 4.8 and 8.4 GHz in the extragalactic radio 
sources PKS 0405-385 (Kedziora-Chudczer et al. 
1997), PKS 1257-326 (Bignall et al. 2003) and 
J1819-I-3845 (Dennett-Thorpe & de Bruyn 2000). 
Weak interstellar scintillation has been demon- 
strated as the principal cause of these rapid vari- 
ations in each of these sources, PKS 0405-385 
(Rickett et al. 2002), PKS 1257-326 (Bignall et al. 
2003) and J1819-F3845 (Dennett-Thorpe & de 
Bruyn 2003), so it follows that WISS should 
also be considered as potentially responsible for 
the variability of the lines in Circinus. 

These three "fast" AGN scintillators exhibit 
such rapid variability because they undergo weak 
scattering and have been found to be behind 
nearby scattering screens. The evidence for weak 
scattering is well established; the flux density vari- 
ations are very broad band and are highly corre- 
lated at 5 and 8 GHz, and each of these sources 
shows much slower variations at 1.4 and 2.4 GHz 



in the strong scattering regime. 

Another strong similarity between the well es- 
tablished AGN weak scintillation and the rapid 
variations in the Circinus masers is the presence of 
a deep first minimum in the ACF. This is clearly 
illustrated in Figure 4, which is the ACF of the 
strongest line in the DOY 124 observations. In 
the case of PKS 0405-385, the deep first mini- 
mum has been shown to indicate scattering in an 
anisotropic ISM (Rickett et al. 2002). 

The evidence for nearby screens for the fast 
scintillators is also compelling. The linear size of 
the first Fresnel zone is proportional to (D/i/)*^-^, 
and hence the corresponding time-scale for travers- 
ing this zone is proportional to (D/z/)"-^ ^^jgj^, 
where D is the screen distance, Vism is the speed 
of the interstellar medium (ISM) and v the fre- 
quency. Thus for the measured values of Vism for 
these sources, ^ 30 km/s, and for the observed 
frequencies of '--^ 5 - 8 GHz where weak scattering 
predominates in AGN, screen distances of ~ 15 
to 30 pc are impUed for PKS 0405-385 (Rickett 
et al. 2002), PKS 1257-326 (Bignall et al. 2003) 
and J1819-h3845 (Dennett-Thorpe & de Bruyn 
2002). 

For Circinus, and following Walker (1998); 

Tehar = 1-2 X lO^X (Vism)-^X {B/vf-\ 

where Tchar is the characteristic time-scale of 
the observed quasi-sinusoidal variations (in sec- 
onds). Substituting Tchar = 700 seconds, typ- 
ical of the time-scales in Table 1, Vism = 50 
kms~^(Rickett et al., 1995), yields a screen dis- 
tance, D, of 20 pc. If the screen is moving at 
the ^ 30 kms~^of the local standard of rest (LSR) 
and the values determined for the fast scintillators, 
then the screen distance is reduced to 7 pc. These 
screen distance estimates are remarkably similar 
to those found for the three fast AGN scintilla- 
tors. At 22 GHz, the angular size of the first 
Fresnel zone is 12 /las for a screen distance of 20 
pc. This is much larger than the expected size of 
the individual maser features. Galactic masers ob- 
served using VLBI subtend milli-arcsecond sizes at 
kpc distances which implies micro-arcsecond sizes 
at Mpc distances as the physical sizes are likely 
to be comparable. The Tidbinbilla observations 
show that all of the lines that have a sufficient 
SNR appear to exhibit rapid variations, as might 
be expected from the large angular size of the first 
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Frcsncl zone. This, in turn, allows an estimate of 
the transverse extent of the scattering screen. The 
overall angular scale of the maser emission from 
the galaxy is ^ 80 mas (Grecnhill et al., 2003), 
which, at 20 pc projects to a linear size of 2.4 
X 10^ km. If all of the disc and outflow lines scin- 
tillate rapidly, then this represents a lower limit to 
the size of the nearby scattering screen. Circinus 
is the only source which has such close multiple 
scintillating components, so this is the first time 
that such a size estimate can be derived. 

For this to be weak scattering requires a tran- 
sition frequency between and weak scattering (uq) 
of somewhat less than 22 GHz. The mean mod- 
ulation index is 0.19, so that following Walker's 
(1998) equation 6, and assuming quite reasonably 
that the maser features are unresolved on a 12 /xas 
scale, yields a transition frequency of 7 GHz, little 
higher than the 3 to 5 GHz found for the AGN, 
albeit at higher Galactic latitudes. We are really 
only looking here at the scattering that is tak- 
ing place at the nearby screen. Circinus lies at a 
low Galactic latitude of -3.8°and so the radiation 
is passing through a dense and turbulent region 
of the Galaxy. Models of the Galactic electron 
distribution such as Cordes and Lazio's NE2001 
model^ (Cordes & Lazio 2003) predicts a transi- 
tion frequency of 56 GHz for the line of sight to 
Circinus, placing the 22 GHz water maser radia- 
tion in the strong scattering regime. Also, as this 
takes the entire path length through the Galac- 
tic plane into account, the screens are likely to 
be at a distance of hundreds or thousands of par- 
sees. We will address the possibility of diffr active 
scintillation in the next section, but these distant 
screens are likely to also cause refractive scintil- 
lation. This will have a characteristic timescale 
of several hours and a relatively large modulation 
index, given by = ^^3, where is the scatter- 
ing strength (Walker 1998). For the predicted 
transition frequency of 56 GHz, this results in a 
modulation index of ~ 0.6. 

There is a clear observational bias present in 
our data towards detecting rapid variations on a 
time scale of minutes to hours but missing those 
variations with time-scales of many hours. Both 
our data and the well-sampled part of Greenhill 
et al.'s data span no more than 3 to 4 hours. 

^Available at http://rsd-wwvj.nrl.navy.mil/7213/lazio/ne-model/ 



However, both the 1995 Parkes and the 1996 Tid- 
binbilla data show evidence of longer-term varia- 
tions with time-scales of at least several hours. It 
would appear highly likely, if WISS is the mecha- 
nism responsible for these variations, longer dura- 
tion observations will uncover longer time-scales, 
of hours to days, which would be caused by scat- 
tering at a more distant screen (s). This would be 
the first time that scattering from multiple screens 
has been observed. 

Using a local screen, we can explain a many of 
the observed characteristics of the variability in 
the Circinus megamasers. However, one problem 
is that we do not see all the lines varying on a 
single timescale which is expected for weak scat- 
tering. The large size of the first Fresnel zone 
means that source structure is unlikely to be re- 
sponsible for the range of observed timescales un- 
less both the source and scintillation pattern are 
highly anisotropic. Axial ratios of 6 or more have 
been suggested for scintillating sources such as 
PKS 0405 385 and J1819-h3845 and VLBI images 
of Galactic masers often show elongated structures 
(Imai, Deguchi & Sasao 2002). Another pos- 
sibility is that refractive scattering at some dis- 
tant screen has broadened the image of the maser 
sources to a size comparable to the first Fresnel 
zone of the local scattering screen. Anistropic 
structures will create distinctive patterns in the 
annual cycle. 

It would appear that the rapid variability ob- 
served in the maser lines in the Circinus galaxy 
is indeed consistent with weak interstellar scintil- 
lation, based on the striking similarities with the 
variations observed in the other three fast AGN 
scintillators. Like the fast AGN scintillators, rapid 
variations imply the presence of a nearby scat- 
tering screen, and the presence of such a nearby 
screen is now more readily acceptable given those 
already found. These screens for the three AGN 
and Circinus are spread widely across the sky and 
suggest the possible presence not so much of a " lo- 
cal bubble" (Bhat, Gupta & Rao, 1998) as the 
remains of a burst bubble. 

3.2. Diffractive scintillation 

Because of the rapidity of the intensity varia- 
tions in Circinus, Greenhill et al. (1997) and Mal- 
oney (2002) both suggested diffractive interstellar 
scintillation (DISS) in the turbulent ISM, rather 
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than the weak scattering in a very nearby screen 
that we have considered above. With DISS it 
is feasible to produce such short time scales for 
the maser lines at 22 GHz, using a distant screen. 
The diffractive timescale is given by tdis = tp/C 
and for a scattering strength of 5 (as predicted by 
the NE2001 model (Cordcs & Lazio 2003)) and 
a screen at 1 kpc, the corresponding diffractive 
timescale is 1000 s. 

However, there are several problems in apply- 
ing diffractive scintillation to the Circinus mega- 
masers. Firstly, as noted by Greenhill et al. and 
Maloney, the observed ~ 700 second time-scale ob- 
served at 22 GHz implied a time scale that was 
significantly shorter than found in pulsars at com- 
parable Galactic latitudes at lower frequencies. 
The diffractive timescale varies with frequency as 
— ^'^ which means that the diffractive timescale at 
1 GHz is a factor of 40 less than at 22 GHz (Walker 
1998). The predicted timescale is therefore on the 
order of 20 seconds, which is comparable to some 
of the most rapidly scintillating pulsars (Cordes 
1986). 

However, Circinus is extragalactic while most 
pulsars are within the Galaxy and as such ex- 
perience less scattering. Also, the combination 
of steep pulsar spectra, narrow diffractive band- 
widths and intrinsic pulse variations makes it ex- 
tremely difficult or impossible to observe very 
short timescale intensity variability in most pul- 
sars. Intrinsic pulse variations have very large 
modulation indices (Kramer et al.2003 reports a 
modulation index of 1.21 for intrinsic variations 
at 1.4 GHz, in the pulsar B1133-I-16) and the vari- 
ations in the pulses are uncorrclatcd. Looking 
for ISS involves averaging over a large number of 
pulses and then examining the variability of the 
averaged data. The averaging time has to be in 
excess of a hundred pulses to remove the effect of 
this intrinsic variability and hence timescales of 
the order of 10 seconds or less cannot be easily 
measured directly. Estimating the strength of the 
scattering could be attempted by measuring dy- 
namic spectra of the two high dispersion measure 
(DM) pulsars (see below) close to the line of sight 
to the Circinus galaxy. 

In modeling the variability of the Circinus 
megamasers as due to diffractive scintillation, we 
note that the masers have a modulation index of 

0.20, rather than unity as expected for diffrac- 



tive scintillation. This implies that the source 
has structure on scales greater than the diffractive 
limit. This extended structure acts as a spatial 
filter, reducing both the speed and amplitude of 
the variations. 

In quenched diffractive scintillation, the mod- 
ulation index is reduced to /x = Os/Odm and the 
timescale is increased by a factor of / ^diff so that 
the dominant lengthscale is now the linear size of 
the source, as projected onto the scattering screen. 
Using this model the observed timescale of a maser 
feature would be proportional to its angular size. 
Taking the most rapidly varying feature with its 
700 seconds timescale and modulation index of 
0.20, we can apply plausible estimates of source 
size and ISM velocity to investigate the screen. 
Using an ISM speed of 50 kms~^and a source size 
of 1 /xas, we obtain a screen distance of ~ 230 pc, 
with a scattering strength of 18 and a transition 
frequency of 120 GHz. Decreasing the assumed an- 
gular size of the masers increases both the implied 
screen distance and scattering strength. For a 0.25 
/xas maser (which corresponds to a linear size of ~ 
1 AU at Circinus), the implied screen distance is 
~ 1 kpc with a scattering strength of ^ 40 and a 
transition frequency of ~ 200 GHz. While this is 
extremely high and well above that predicted by 
the NE2001 model for the line of sight to Circinus, 
it is characteristic of lines of sight at lower Galac- 
tic latitudes. We have inspected the ATNF pulsar 
database^ and it is worth noting that the two pul- 
sars closest to Circinus both have quite high DMs 
~ 250 pc/cm"^, which is higher than that of other 
pulsars at comparable Galactic latitudes. We have 
used the ATNF pulsar database to produce a DM 
map of the Galactic plane near to the line of sight 
to Circinus (Figure 5), which shows that the line 
of sight to Circinus is associated with high DMs. 
This suggests that either the pulsars in this re- 
gion arc more distant than some of the other, less 
scattered pulsars or that there is some enhanced 
scattering occuring along that line of sight, or per- 
haps both. At present, none of the pulsars used in 
the DM map have had their distances measured in- 
dependently. Ha images of the same region show 
strong emission nearby, which is suggestive that 
enhanced scattering may be present. 

Another problem with applying a diffractive 

^ http://www. atnf. csiro. au/research/pulsar/psrcat/ 
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scintillation model is that the expected anticor- 
relation between timescale and modulation index 
is not apparent in the data. Figure 3 illustrates 
this well with the very large amplitude flux varia- 
tion in the 295 kms~^line occurring on a very long 
timescale, contrasted with the rapid, lower ampli- 
tude variations in the 256 kms~"'^linc. This may 
be an effect of source and/or screen anisotropy 
which is suggested by the deep first minima of the 
ACFs, which will produce a range of timcscales 
depending on the relative orientations of the axes 
of anisotropy and screen velocity. 

4. Conclusions 

We have shown that the variability of the Circi- 
nus megamasers has continued since its first detec- 
tion and that it can be explained through ISS. The 
Circinus megamasers offer a unique opportunity to 
study both the interstellar scattering medium and 
the structure of extragalactic megamasers at high 
angular resolution. 

There are several experiments that can be made 
to confirm that scintillation is the cause of the 
variability and identify whether is due to a weakly 
scattering local screen or a strong distant one. 

Monitoring the variability of the Circinus mega- 
masers will allow us to test for the presence of an 
annual cycle in its characteristic timescale. The 
annual cycle is created by the motion of the Earth 
around the Sun and can be used to measure the pe- 
culiar velocity of the scattering material, estimate 
the extent and orientation of anisotropics in the 
medium and to create a model of the source itself 
as has been accomplished for PKS 0405-385 (Rick- 
ett et al. 2002), J1819-F3845 (Dennett-Thorpe & 
de Bruyn 2003) and PKS 1257 326 (Bignall et al. 
2004). Both of our ISS scenarios expect that 
the scattering medium will be anisotropic and 
that source structure is directly involved in the 
timescale of variability. This means that different 
maser features may display different annual cycles 
as a result of their structure. In the case of a local 
screen, its intrinsic velocity is likely to be compa- 
rable to local standard of rest, producing a clear 
annual cycle with speeds < 50 kms~^. A more 
distant screen may have some significant intrinsic 
velocity which will make the annual cycle harder to 
detect as the difference between the high and low 
speeds becomes proportionally less. The relative 



transverse speed of the ISM can also be directly 
measured by observing the time delay between the 
scintillation patterns recorded at two widely sep- 
arated telescopes as has been done for PKS 0405- 
385 (Jauncey et al. 2000), PKS 1257-326 (Bignall 
et al. 2004) and J1819-I-3845 (Dennett-Thorpe & 
de Bruyn 2002), using intercontinental baselines. 
A direct measurement of the ISM velocity is an 
extremely valuable constraint in producing mod- 
els of anisotropic structure. 

These experiments will not conclusively distin- 
guish between a local or distant screen model. One 
test would be to investigate longer timescale vari- 
ability, which is predicted by both models as a 
result of refractive scintillation at some distant 
screen. In the local screen scenario, the refrac- 
tive timescale would be decoupled from that of the 
rapid variations and while in the distant screen, 
both the diffractive and refractive variations would 
display the same annual cycle. The maser s should 
all be point sources to refractive scintillations and 
so the timescale and modulation index is expected 
to be the same for all the masers. Also, the pul- 
sars near to Circinus provide another probe of 
the ISM. If characteristic timescales and diffrac- 
tive bandwidths can be measured from dynamic 
spectra then these can be used to make estimates 
of the screen strength and distance. 

We have shown that the presence of variability 
in the Circinus megamasers can be explained as 
due to ISS and have presented two models. At 
this stage, both require unusual, although not un- 
precedented, conditions to explain some charac- 
teristics of the Circinus megamasers. The local 
screen model using weak scintillation can explain 
the observed variability using a local, weakly scat- 
tering screen with a distance < 20 pc and i>q ~ 7. 
This is comparable to what has been observed in 
the most rapidly scintillating AGN sources. The 
appearance of the light curves is also strongly rem- 
iniscent of those seen in scintillating AGN. How- 
ever, it does have some problems in explaining 
the observed range of timescales, as the angular 
size of the first Fresnel zone vastly exceeds that 
of any realistic source, unless there is also scat- 
tering by a more distant screen. Also, the low 
Galactic latitude of Circinus and the predictions 
of models of the Galactic electron distribution sug- 
gest that strong scattering and distant screens 
should be involved. The distant, strongly scatter- 
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ing screen model we have proposed satisfies most 
of our observations and fits with the predictions 
of electron density models. However, it predicts 
an anti-correlation between modulation index and 
timescale which we do not observe. Further obser- 
vations will allow us to refine our models and to 
resolve some of our present uncertainties. 

Acknowledgements 

The ATNF is funded by the Commonwealth 

Government for operation as a National Facility 
by CSIRO. Financial support for this work was 
provided by the Australian Research Council. We 
also wish to thank Edward King for making the 
Tidbinbilla observations and developing the soft- 
ware to read in the data, and to our referee for 
their extremely helpful comments. 

REFERENCES 

Bhat, N D R., Gupta, Y., Rao, A. P, 1998, ApJ, 
500, 262 

Bignall, H. E., Jaunccy, D. L., Lovell, J. E. J., 
Tzioumis, A. K., Kedziora-Chudczer, L., Mac- 
quart, J. -P., Tingay, S. J., Rayner, D. P., Clay, 
R. W., 2003, ApJ, 585, 653 

Bignall, H. E., Jauncey, D. L., Lovell, J. E. J., 
Macquart, J. -P., Tzioumis, A. K., Rickctt, B. 
J., Ojha, R., Carter, S., Cimo, G., EUingsen, S., 
McCulloch, P. M., 2004, Proceedings of the 7*'' 
EVN Symposium, Bachiller, Colomer, Desmur 
fc Vicente Eds 

Cordcs, J. M., 1986, ApJ, 311, 183 

Cordes, J. M., Lazio, T. J. W., 2003, arXiv 
preprint - arXiv:astro-ph/0207156 v3 

Dennett-Thorpe, J., de Bruyn, A. G., 2000, ApJ, 
529, L65 

Dennett-Thorpe, J., de Bruyn, A. G., 2002, Na- 
ture, 415, 57 

Dennett-Thorpe, J., de Bruyn, A. G., 2003, A&A, 

404, 113 

Edelson, R. A., Krolik, J. H., 1988, ApJ, 333, 646 

Gardner, F. F., Whiteoak, J. B. 1982, MNRAS, 
201 P13 



Greenhill, L. J., EUingsen, S. P., Norris, R. P., 
Gough, R. G., Sinclair, . W., Moran, J. M., 
Mushotzsky, R., 1997, ApJ, 474, L103 

Greenhill, L. J., Booth, R. S., EUingsen, S. P., 
Herrnstein, J. R., Jauncey, D. L., McCulloch, 
P. M., Moran, J. M., Norris, R. P., Reynolds, 
J. E., Tzioumis, A. K., 2003, ApJ, 590, 162 

Imai, H., Deguchi, S., Sasao, T., 2002, ApJ, 567, 
971 

Jaimcey, D. L., Kedziora-Chudczer, L. L., Lovell, 
J. E. J., Nicolson, G. D., Perley, R. A., 
Reynolds, J. E., Tzioumis, A. K., Wieringa, M. 
H., 2000, in Proc. lAU 97, Extragalactic Ra- 
dio Sources, ed. D. S. Heeschen &: C. M. Wade, 
(Dordrecht: Kluwer), 327 

Jauncey, D. L., Bignall, H. E., Lovell, J. E. J., 
Kedziora-Chudczer, L., Tzioumis, A. K., Mac- 
quart, J.-P., Rickett, B. J., 2003, in ASP Conf. 
Ser. 300, Radio Astronomy at the Fringe, ed. J. 
A. Zensus, M. H. Cohen, & E. Ros (San Fran- 
cisco: ASP), 199 

Kedziora-Chudczer, L., Jauncey, D. L.. Wieringa, 
M. H., Walker, M. A., Nicolson, G. D., 
Reynolds, J. E., Tzioumis, A. K., 1997, ApJ, 
490, L9 

Kramer, M., Karastergiou, A., Gupta, Y., John- 
ston, S., Bhat, N. D. R., Lyne, A. G., 2003, A 
& A, 407, 655 

Macquart, J.-R, Jauncey, D. L., 2002, ApJ, 572, 
786 

Maloney, P.,R., 2002, PASA, 19, 401 

Rickett, B. J., 1990, Annual Review of Astronomy 

and Astrophysics, 28, 561 

Rickctt, B. J., 2002, PASA, 19, 100 

Rickett, B. J., Kedziora-Chudczer, L., Jauncey, D. 
L., 2002, ApJ, 581, 103 

Rickett, B. J., Quirrenbach, A., Wegner, R., 
Krichbaum, T. P., Witzel, A., A&A 1995, 293, 
479 

Walker, M. A., 1998, MNRAS, 294, 307 



This 2-column preprint was prepared witii the AAS lA'H^X 
macros v5.2. 



7 



3.5 



Systemic Velodty 439 kms" 



Heliocentric Velocity {kms~ ) 



Fig. 1 . — Average spectrum of the Circinus mega- 
masers, made using the Tidbinbilla telescope. The 
two observations made on DOYs 124 and 148, 
of the red- and blue-shifted emission respectively, 
have been amalgamated. The velocity range cov- 
ered by both observations is between 402 and 490 
kms~^. 
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Fig. 3. — Light curves of the two strongest lines 
observed on day 148 of 1996 from Tidbinbilla. The 
upper panel shows the light curve of the 256.79 
kms~^line and the lower panel the light curve of 
the 295 kms~^line. The two lines display very 
different timescales of variability. While the 256 
kms~^line completes two quasi-periods, the 295 
kms~^line shows perhaps half of one. 
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Fig. 2. — Light curves of the three strongest lines, 
day 124 of 1996 from Tidbinbilla. A range of 
timescales is apparent in the data, with short- 
timescale, low amplitude variations overlaid on a 
much larger, slower trend in the 558 kms~'^line. 
The strongest line, at 564 kms~^, is observed to 
more than halve its flux density in less than 10 
minutes. This is the largest proportional change 
observed in this set of observations. 
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Fig. 4. — Autocorrelation function of the 564 
kms~^line. This was computed using the discrete 
autocorrelation function as described in Edelson & 
Krolik (1988). The lags are binned in 40 second 
intervals. Data are indicated by dots and verical 
bars indicate the one a estimation error, following 
Rickett et al. (2002). The large estimation errors 
are due to the short observing interval. Missing 
lags correspond to where the sampling function 
has left too few data to make a reliable estimation 
of the ACF. The ACF shows a deep first minimum 
and a strong quasi-periodic character. 
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Fig. 5. — This is a smoothed map of DM for 
pulsars along the Galactic plane near Circinus. 
The position of the Circinus Galaxy is indicated 
with a white star while the individual pulsars 
are indicated with black ones. The data was 
smoothed using two-dimensional spline interpola- 
tion. The data has been taken from the ATNF 
pulsar database 
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Table 1 

\/^4riability parameters of the obseraable lines 



Velocity 


FWHM 


< S > 


(T 




Tchar 


First min. 


(kms-^) 


(kms-^) 


(Jy) 


(Jy) 




(Hours) 


Depth 




Parkes, day 289 of 1995 


295.03 ± 0.02 


0.59 ± 0.04 


11.97 


1.22 


0.1 


0.07 ± 0.01 


-0.4 ± 0.2 


556.97 ± 0.05 


1.32 ± 0.16 


11.36 


1.51 


0.13 


0.4 ± 0.2 


No min. 


562.05 ± 0.03 


1.04 ± 0.08 


18.03 


3.64 


0.20 


0.3 ± 0.2 


-0.2 ± 0.3 


Tidbinbilla, day 124 of 1996 


514.74 ± 0.04 


1.28 ± 0.08 


2.41 


0.34 


0.42 


0.4 ± 0.2 


-0.4 ± 0.4 


557.03 ± 0.03 


0.86 ± 0.06 


4.93 


0.75 


0.15 


0.3 ± 0.1 


-0.5 ± 0.4 


558.13 ± 0.02 


0.87 ± 0.04 


7.66 


1.14 


0.15 


0.5 ± 0.3 


No min. 


560.3 ± 0.3" 


1.9 ± 0.4 


4.1 


0.4 








563.24 ± 0.05 


1.88 ± 0.07 


5.67 


0.69 


0.12 


0.26 ± 0.13 


-0.8 ± 0.3 


564.22 ± 0.05 


0.84 ± 0.05 


11.32 


1.84 


0.16 


0.09 ± 0.02 


-0.4 ± 0.2 


578.46 ± 0.03 


0.78 ± 0.11 


1.59 


0.37 


0.23 


0.16 ± 0.06 


-0.8 ± 0.3 


598.17 ± 0.03 


1.18 ± 0.13 


2.26 


0.23 


0.10 


0.16 ± 0.06 


-0.2 ± 0.3 


614.53 ± 0.03 


0.82 ± 0.08 


1.96 


0.24 


0.12 


0.09 ± 0.02 


-0.8 ± 0.2 


642.65 ± 0.09 


0.96 ± 0.18 


0.91 


0.25 


0.27 


0.24 ± 0.1 


-0.2 ± 0.3 


Tidbinbilla, day 148 of 1996 


256.79 ± 0.02 


0.79 ± 0.10 


2.50 


0.34 


0.13 


0.10 ± 0.03 


-0.6 ± 0.25 


295.00 ± 0.01 


0.63 ± 0.02 


5.68 


1.18 


0.21 


0.34 ± 0.2 


No min. 



*This is a blended feature with two broad components which cannot be resolved. 
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